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LONG-TERM  GOALS 

Advancing  aerosol  optical  property  studies  from  both  observational  and  modeling  perspectives  is 
necessary  for  military  applications,  especially  for  operations  with  visibility-sensitive  instruments  such 
as  electro-optical  (EO)  systems.  The  long-term  goal  of  this  study  is  to  fully  evaluate  existing  space- 
borne  aerosol  optical  property  data  sets,  innovatively  construct  new  satellite  aerosol  retrieval  schemes, 
and  develop  advanced  methods  for  assimilating  existing  and  new  multi-sensor  aerosol  measurements 
for  improving  Navy’s  visibility  and  EO  propagation  forecasting  capabilities. 

OBJECTIVES  (abstract  from  proposal) 

Critical  to  both  military  and  civilian  applications,  the  Navy  Aerosol  Analysis  and  Prediction  System 
(NAAPS)  is  the  world's  only  truly  operational  global  aerosol  and  visibility  forecasting  model.  Recent 
studies  indicate  that  the  assimilation  of  satellite  observations  significantly  improves  NAAPS  aerosol 
forecasting  capability  and  reliability.  Io  fully  utilize  the  wide  breadth  and  depth  of  various  current 
satellite  observations  and  prepare  for  future  reductions  in  aerosol  sensing  satellites  over  the  next 
decade,  we  propose  to  construct  a  multi-channel,  multi-sensor,  and  multi-task  assimilation  system  to 
improve  NAAPS  forecasts  for  both  current  and  future  applications.  The  specific  objectives  are  to: 

1.  Finalize  over-land  and  over-ocean  aerosol  assimilation  methods  using  operational  data 
assimilation  (DA)  quality  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  and 
Multiangle  Imaging  SpectroRadiometer  (MISR)  aerosol  products,  and  develop  a  framework  for 
considering  current  and  future  satellite  aerosol  products. 

Develop  forward  models  to  enable  a  radiance  assimilation  capability  by:  1)  improving  forecast 
performance  over  cloudy  regions  using  the  Ozone  Monitoring  Instrument  (OMI)  Aerosol 
Index;  and  2)  preparing  for  the  post-MODIS/MISR  era  using  the  Geostationary  Operational 
Environmental  Satellite  (GOES). 

3.  Improve  model  representation  of  aerosol  vertical  profiles  and  the  accuracy  of  aerosol  speciation 
in  NAAPS  through  the  use  of  a  3-D  aerosol  assimilation  method  on  satellite-based  lidar 
observations  and  a  generalized  Angstrom  exponent  assimilation  scheme. 

4.  Develop  an  improved  3-D  parameterization  for  satellite  observation  and  model  forecasting 
error  matrices  using  ground  observations  from  the  NASA  Aerosol  Robotic  Network 
(AERONET)  and  the  Micropulse  Lidar  Network  (MPLNET). 
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APPROACH 


To  achieve  the  proposed  objectives  and  long-term  goals,  research  approaches  of  the  project  include 
constructing  DA-grade  satellite  data  for  aerosol  DA,  developing  new  and  advanced  aerosol 
assimilation  schemes,  and  advancing  new  and  cutting-edge  research  to  assist  the  proposed  research 
topics. 

Supported  by  this  grant,  Shi  et  al.  (2011)  suggested  that  large  spatial  discrepancies  exist  in  operational 
satellite  aerosol  products  such  as  MODIS  and  MISR.  Thus,  before  applying  satellite  aerosol  retrievals 
in  aerosol  analysis  and  forecasts,  the  uncertainties  and  biases  in  operational  satellite  aerosol  products 
need  to  be  fully  evaluated.  In  the  past  few  years  of  the  project  period,  MODIS  Deep  Blue  (DB)  and 
MISR  aerosol  products  have  been  studied,  and  schemes  for  constructing  DA-grade  MODIS  DB  and 
MISR  aerosol  products  have  been  developed  and  transitioned  to  the  Naval  Research  Laboratory 
(NRL).  Continuing  from  FY14’s  effort,  in  FY15,  the  newly  released  Collection  6  MODIS  Dark  Target 
(I)T)  over  water  products  have  been  evaluated  for  their  potential  applications  in  operational  aerosol 
DA  (Shi  et  al.,  2014). 

New  improvements  are  also  made  in  developing  advanced  aerosol  DA  schemes.  In  the  past  few  years 
of  the  study  period,  a  3-D  CALIOP  DA  scheme,  which  is  coupled  with  2D  aerosol  optical  depth 
(AOD)  assimilation  using  MODIS  and  MISR  aerosol  data,  has  been  completed  (Zhang  et  al.  2011; 
2014)  and  transitioned  to  NRL.  In  F  Y14,  a  prototype  GOES  radiance  assimilation  system,  specifically 
lor  aerosol  analysis  and  forecasts,  was  designed  and  developed  through  coupling  of  the  Community 
Radiative  Transfer  Model  (CRTM)  with  NAAPS  and  NAVDAS.  In  FY15,  attempts  have  been  made 
to  assimilate  real  GOES-13  data  (e.g..  Fig.  1). 

Lastly,  in  FY15,  exploratory  studies  have  continued  in  research  directions  that  could  benefit  the  overall 
research  goal  of  the  project.  A  revised  nighttime  aerosol  retrieval  scheme  (based  on  Johnson  et  al., 
2013,  developed  with  the  support  of  this  project)  has  been  developed  by  using  the  variance  in 
radiances  within  an  artificial  light  source  to  retrieve  AOD  (McHardy  et  al.,  2015).  Research  efforts 
have  also  continued  to  study  the  spatial  and  temporal  distributions  of  above  cloud  aerosol  events 
(Alfaro-Contreras  et  al.,  2015),  as  well  as  complete  the  analysis  of  the  effect  of  sub-surface  bubbles  on 
aerosol  retrievals  (Christensen  et  al.,  2015),  which  has  not  been  considered  in  aerosol  retrievals  prior  to 
our  study. 


WORK  COMPLETED 

Overall  research  progress  of  the  project 

for  the  past  few  years,  this  project  has  focused  on  developing  advanced  aerosol  DA  schemes  and 
constructing  reliable  data  streams  for  operational  aerosol  DA.  Data  source-wise,  DA-grade  Moderate 
Resolution  Imaging  Spectroradiometer  (MODIS)  DeepBlue  (DB)  and  Multiangle  Imaging 
SpectroRadiometer  (MISR)  aerosol  products  have  been  developed.  The  enhanced  southern  ocean 
anomaly  (ESOA),  an  anomaly  found  in  the  MODIS  and  MISR  aerosol  products,  has  been  studied  with 
respect  to  cloud  contamination  and  subsurface  oceanic  bubbles.  Methods  have  also  been  developed  for 
cloud-clearing  of  MISR  observations  with  MODIS  data.  DA  scheme-wise,  continuing  from  the  Pi's 
ONR  YIP  project,  a  2-D/3-D  aerosol  assimilation  package  that  assimilates  both  column-integrated 
aerosol  optical  depth  (AOD)  data  and  lidar  vertical  aerosol  extinction  profiles  from  Cloud-Aerosol 
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Lidar  and  Infrared  Pathfinder  Satellite  Observations  (CALIOP)  has  been  completed  and  transitioned  to 
NRL.  A  new  prototype  model  for  directly  assimilating  GOES  radiance  has  also  been  developed.  The 
regional  and  global  impact  of  assimilating  MODIS  Dark  Target  (DT),  MODIS  DB,  MISR  and 
CALIOP  aerosol  products  for  aerosol  analysis  and  forecasts  has  been  evaluated.  Beyond  the  proposed 
research  topics,  new  and  cutting-edge  research  efforts  that  directly  relate  to  the  proposed  theme  have 
been  explored.  In  particular,  innovative  methods  have  been  developed  for  retrieving  nighttime  AOD 
from  Visible  Infrared  Imaging  Radiometer  Suite  (VIIRS)  observed  artificial  light  sources.  The  spatial 
and  temporal  variations  of  above  cloud  aerosol  events  are  evaluated  and  the  impacts  of  above  cloud 
aerosols  to  cloud  property  retrievals  are  studied.  In  addition,  the  relationship  between  surface 
particulate  matter  concentration  and  column  integrated  aerosol  optical  thickness  has  been  investigated 
as  functions  of  data  quality  and  aerosol  vertical  distribution.  Research  findings  from  this  project  have 
been  reported  through  a  total  of  more  than  25  submitted/published  authored  and  coauthored  peer- 
review  journal  papers  and  more  than  30  conference  presentations. 

Research  progress  for  the  last  reporting  period  (Oct.  1  201 4  -Mav  30.  201  5) 

Continuing  from  FY14’s  effort,  a  revised  version  of  our  nighttime  aerosol  retrieval  method  has  been 
developed.  Currently,  no  reliable  nighttime  aerosol  retrievals  from  passive-based  sensors  are  available 
for  operational  aerosol  forecasts  (Zhang  et  al.,  2011).  While  CALIOP  has  the  capability  of  detecting 
nighttime  aerosol  vertical  distribution,  the  profiling  swath  is  limited.  In  FY13,  a  spatial  contrast 
method  was  developed  by  deriving  AOD  through  evaluating  the  radiance  contrast  between  an  artificial 
light  source  and  nearby  background  pixels  (Johnson  et  al.,  2013).  A  revised  method  has  also  been 
attempted  by  retrieving  nighttime  aerosol  optical  thickness  through  measuring  the  variance  in  radiance 
within  an  artificial  light  source  (McHardy  et  al.,  2015).  The  revised  method,  which  excludes  the  need 
for  nearby  pixels,  has  the  potential  to  be  applied  to  a  much  larger  domain.  Validated  against  nighttime 
High  Spectral  Resolution  Lidar  (HSRL)  measurements  and  adjacent  daytime  AERONET  data,  the 
retrievals  from  the  proposed  methods  are  very  promising  (McHardy  et  al.,  2015). 

In  FY14,  a  GOES  radiance  assimilation  package  was  developed  and  tested  with  synthetic  data. 
Implementation  of  the  GOES  radiance  assimilation  package  with  the  use  of  real  GOES- 13  data  was 
attempted  in  FY15.  A  total  of  3  months  of  GOES- 13  data  from  June-Aug.  2007  were  processed  and 
cloud  screened.  Collocated  GOES-13  and  AERONET  data  are  used  to  evaluate  the  simulated  GOES 
radiance  from  a  GOES  radiance  simulator  (a  fully  coupled  NAAPS  and  CRTM  system).  The  gridded 
GOES  radiance  data  are  further  used  for  testing  the  GOES  radiance  assimilation  system  (Fig.  1). 

In  early  2014,  Collection  6  Aqua  MODIS  aerosol  products  were  released.  We  have  evaluated  over¬ 
water  Collection  6  MODIS  D  T  products  with  respect  to  cloud  contamination,  aerosol  microphysical 
bias,  and  near  surface  wind  conditions  (Shi  et  al.,  2014).  The  full  evaluation  of  the  DB  and  DT  Aqua 
and  Terra  MODIS  products  will  likely  be  a  future  study  as  the  Collection  6  Terra  MODIS  products 
have  been  released  one  month  (April  2015)  prior  to  the  end  of  the  project. 

Lastly,  attempts  have  also  been  made  to  study  the  spatial  and  temporal  variations  in  aerosol  vertical 
distribution.  While  we  are  still  working  on  converting  the  study  to  a  peer-reviewed  journal  paper,  the 
first  phase  of  the  study  has  been  presented  at  the  2014  Fall  AGU  meeting  (Toth  et  al.  2014). 
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RESULTS 


An  improved  method  for  retrieving  nighttime  aerosol  optical  thickness  from  the  VIIRS  Dav/Nipht. 
Band  (McHardv  et  al.  2015.  submitted) 

Using  Visible/Infrared  Imaging  Radiometer  Suite  (VIIRS)  Day/Night  Band  (DNB)  data,  a  method, 
dubbed  the  “variance  method"  was  developed  for  retrieving  nighttime  aerosol  optical  thickness  (x) 
values  through  the  examination  of  the  dispersion  of  radiance  values  above  an  artificial  light  source. 
Based  on  improvement  to  a  previous  algorithm,  this  updated  method  derives  a  semi-quantitative 
indicator  of  nighttime  x  using  artificial  light  sources.  Nighttime  x  retrievals  from  the  newly-developed 
method  are  inter-compared  with  an  interpolated  value  from  late  afternoon  and  early  morning  ground 
observations  from  four  AErosol  RObotic  NETwork  (AERONET)  sites  as  well  as  column-integrated  x 
from  one  High  Spectral  Resolution  Lidar  (HSRL)  site  at  Huntsville,  AL  during  the  NASA  Studies  of 
Emissions  and  Atmospheric  Composition,  Clouds  and  Climate  Coupling  by  Regional  Surveys 
(SEAC  RS)  campaign,  providing  diurnal  coverage.  Sensitivity  studies  were  performed  to  examine  the 
effects  of  lunar  illumination  on  VIIRS  x  retrievals  made  via  the  variance  method,  revealing  that  lunar 
contamination  may  have  a  smaller  impact  than  previously  thought.  However,  the  small  sample  size  of 
this  study  limits  the  conclusiveness  thus  far.  VIIRS  x  retrievals  yield  a  coefficient  of  determination  (r2) 
of  0.60  and  a  root-mean-squared-error  (RMSE)  of  0.18  when  compared  against  straddling  daytime- 
averaged  AERONET  x  values.  Preliminary  results  suggest  that  artificial  light  sources  can  be  used  for 
estimating  regional  and  global  nighttime  aerosol  distributions  in  the  future. 

Evaluating  MODIS  collection  6  dark  target  over  water  aerosol  products  for  multi-sensor  data  fusion 

(Shi  et  al.  2014;  presented  at  the  2014  AGU  fall  annual  meeting! 

MODIS  aerosol  products  have  been  widely  used  in  aerosol  related  climate,  visibility,  and  air  quality 
studies  for  more  than  a  decade.  Recently,  the  MODIS  Collection  6  (C6)  aerosol  products  from 
MODIS -Aqua  have  been  released.  The  reported  changes  between  C5  and  C6  include  updates  to  the 
retrieval  algorithms  and  a  new  cloud-filtering  process  for  over-ocean  products.  Thus  it  is  necessary  to 
fully  evaluate  the  C6  products  for  applications  that  require  high  quality  MODIS  aerosol  optical  depth 
data,  such  as  operational  aerosol  data  assimilation.  The  uncertainties  in  the  MODIS  C6  DT  over-ocean 
products  are  studied  through  both  inter-comparing  with  MISR  aerosol  products  and  by  evaluation 
against  ground  truth.  Special  attention  is  given  to  the  low  bias  in  MODIS  DT  products  due  to  the 
mi  ^classifications  of  heavy  aerosol  plumes  as  clouds.  Finally,  a  quality-assured  DA  grade  aerosol 
optical  product  is  constructed  for  aerosol  data  assimilation  related  applications. 

A  Global  and  Regional  Trend  Study  of  the  Vertical  Distribution  of  Aerosols  as  Observed  by  CALIOP 

(Toth  et  al,  2014;  presented  at  the  2014  AGU  fall  annual  meeting) 


Trends  in  atmospheric  aerosol  particle  loading  have  gained  increased  attention  in  recent  years  due  to 
their  impact  on  Earth’s  radiation  budget  and  global  climate  change.  Past  studies  have  examined  this 
topic  through  aerosol  optical  depth  (AOD)  observations  derived  from  passive  satellite  sensors  such  as 
MODIS  and  MISR.  However,  such  passive  sensors  acquire  only  column-integrated  measurements  and 
thus  can  provide  no  insight  into  the  vertical  distribution  of  any  AOD  trends  they  might  detect.  Yet 
knowledge  of  aerosol  vertical  distribution  and  trends  are  critical  for  studies  involving  aerosol  climate 
impacts  and  air  quality.  Using  seven  and  a  half  years  (June  2006  -  December  2013)  of  aerosol  profile 
data  from  Cloud-Aerosol  Lidar  with  Orthogonal  Polarization  (CALIOP),  we  investigate  trends  in 
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CALIOP-derived  AOD  at  various  altitudes  throughout  the  atmospheric  column.  This  differs  from 
existing  CALIOP-based  climatological  studies  that  focus  solely  on  the  mean  state  of  aerosol  vertical 
distribution.  Daytime  and  nighttime  analyses  are  decoupled  to  account  for  differences  in  the  vertical 
structure  of  the  troposphere,  and  detection  sensitivity  of  CALIOP,  between  the  two  regimes.  The 
results  of  this  study  are  presented  globally  and  for  selected  regions.  Of  particular  interest  are  North 
Africa  and  Asia,  as  these  are  areas  with  frequently  high  AOD  and  have  also  exhibited  trends  from 
MODIS  and  MISR  observations. 

Progress  with  GOES  radiance  assimilation 

In  FY15,  the  prototype  system  for  GOES  radiance  assimilation,  which  was  built  in  FY14,  has  been 
tested  with  GOES- 13  data.  Figure  la  shows  the  gridded  (l°xl°  Latitude,  Longitude)  GOES  radiance 
data  from  cloud  free  regions.  Cloud-clearing  of  GOES  data  is  performed  based  on  the  combination  of 
visible/thermal  IR  threshold  methods  and  a  visible  spatial  variance  test.  Figure  lb  shows  the  NAAPS 
aerosol  optical  depth  fields  without  data  assimilation  and  Fig.  lc  shows  the  NAAPS  AOD  fields  after 
radiance  assimilation.  A  biomass-burning  event  in  Central  America,  which  is  not  found  from  the 
NAAPS  natural  run,  is  captured  from  the  GOES  radiance  data  assimilation.  Still,  this  practice  reveals 
potential  problems  in  GOES  radiance  assimilation,  which  includes  cloud  and  glint  contamination  of 
GOES  data,  as  well  as  GOES  visible  channel  signal  degradation.  These  issues  will  be  explored  in  a 
follow  up  ONR  project. 


IMPACT/APPLICATIONS 

A  new  nighttime  aerosol  retrieval  scheme  has  been  developed  with  the  use  of  the  radiance  variance 
within  artificial  light  sources.  The  new  scheme,  which  doesn’t  need  a  priori  knowledge  of  aerosol 
properties  or  background  information,  has  the  potential  to  be  implemented  on  a  regional  domain,  and 
possibly  a  nighttime  aerosol  product  for  aerosol  data  assimilation. 

TRANSITIONS 

None. 
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Figure  1.  a).  Gridded  GOES-13  visible  radiance  (mW/m^sr-cm'1)  over  cloud  free  regions  for  18UTC, 
May  15,  2007.  b)  AOD  (0.55  pm)  distributions  from  NAAPS  natural  run  for  18UTC,  May  15,  2007. 
c).  NAAPS  analysis  with  GOES  radiance  assimilation  using  data  from  Figure  la  (AODs  at  the  0.55 
fim  spectral  channel),  d)  Simulated  (using  the  NAAPS-CRTM  system)  versus  measured  over  water 
GOES- 13  reflectance  using  a  collocated  GOES,  NAAPS  and  AERONET  dataset  for  the  period  of  June 
1  -Aug.  31,  2015.  AERONET  AOD  values  and  NAAPS  aerosol  vertical  distributions  are  used  in  the 
simulations,  e)  VIIRS  retrieved  nighttime  AOD  for  a  new  method  developed  in  this  study  period  as  a 
function  of  straddling  daytime-averaged  AERONET  AOD  for  Huntsville  (square),  Alta  Floresta 
(circle)  and  Grand  Forks  (“x”).  One-to-one  (dotted)  and  best-fit  (solid)  lines  are  also  shown,  f)  VIIRS 
retrieved  AOD  as  a  function  of  Wisconsin  HSRL  retrieved  total-column  AOD  for  Huntsville, 
Alabama. 
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